1. Joint Modeling
1.1 Trace plots for convergence check
The current MCMC setting is:
- 100,000 iteration;
- 90,000 burn-in;
- 10 thinning.
1.2 Gelman and Rubin’s convergence check
## Potential scale reduction factors:
##
## Point est. Upper C.I.
## HDevsum 1.00 1.01
## LDevsum 1.00 1.01
## dh0 1.06 1.16
## dh1 1.05 1.13
## dl0 1.10 1.27
## dl1 1.19 1.50
## dl2 1.02 1.04
##
## Multivariate psrf
##
## 1.2
1.3 ACF Plots
Here we plotted ACF plots for the following variables:
- Total deviance;
- Variables that didn’t pass the convergence check.
1.4 WAIC results
| LevelH | LevelL | |
|---|---|---|
| DIC | 1473.9832 | 22568.0036 |
| DIC3 | 1312.3190 | 22441.7809 |
| PWAIC | 122.9428 | 136.6518 |
| WAIC | 1407.7695 | 22448.2204 |
2. Separate Modeling
2.1 Trace plots for convergence check
The current MCMC setting is:
- 100,000 iteration;
- 90,000 burn-in;
- 10 thinning.
2.2 Gelman and Rubin’s convergence check
## Potential scale reduction factors:
##
## Point est. Upper C.I.
## HDevsum 1.01 1.04
## LDevsum 1.00 1.00
## dh0 2.60 4.97
## dh1 2.57 6.22
## dl0 1.19 1.47
## dl1 1.05 1.14
## dl2 1.03 1.06
##
## Multivariate psrf
##
## 2.3
2.3 ACF Plots
Here we plotted ACF plots for the following variables:
- Total deviance;
- Variables that didn’t pass the convergence check.
2.4 WAIC results
| LevelH | LevelL | |
|---|---|---|
| DIC | 1368.04979 | 30293.5799 |
| DIC3 | 1241.80368 | 30384.2713 |
| PWAIC | 88.31523 | 140.1374 |
| WAIC | 1309.32686 | 30388.2600 |